format short e 
echo on
% Dati da approssimare
xnodi = [ 0 1 2 3 4 5];
fnodi = [10 25 51 66 97 118];
%grafico dei dati
plot(xnodi,fnodi,'*')
hold on
%% costruiamoci la retta di approssimazione nel senso dei minimi 
%% quadrati (retta di best-fitting)
p = polyfit(xnodi,fnodi,1);
xx = linspace(0,5,100);
yp = polyval(p,xx);
pause
% grafico della retta interpolante
plot(xx,yp,'r-')
pause
hold off


% dati relativi alla conducibilita' termica del ferro
xnodi = [ 100 200 300 400 500 600 700 800 900 1000];
fnodi = [1.32 0.94 0.835 0.803 0.694 0.613 0.547 0.487 0.433 0.38];
%grafico dei dati
plot(xnodi,fnodi,'*')
hold on
%% costruiamoci la retta di approssimazione nel senso dei minimi 
%% quadrati (retta di best-fitting)
p1 = polyfit(xnodi,fnodi,1);
xx = linspace(0,1000,100);
yp1 = polyval(p1,xx);
pause
% grafico della retta 
plot(xx,yp1,'r-')
pause
hold on
%% costruiamocial polinomio di 2 grado di approssimazione nel senso dei minimi 
%% quadrati 
p2 = polyfit(xnodi,fnodi,2);
xx = linspace(0,1000,100);
yp2 = polyval(p2,xx);
pause
% grafico 
plot(xx,yp1,'r-',xx,yp2,'g+')
pause
%% costruiamocial polinomio di 3 grado di approssimazione nel senso dei minimi 
%% quadrati 
p3 = polyfit(xnodi,fnodi,3);
xx = linspace(0,1000,100);
yp3 = polyval(p3,xx);
pause
% grafico 
plot(xx,yp1,'r-',xx,yp2,'g+',xx,yp3,'k*')
pause


hold off


